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ABSTRACT 
Ships, submarines, and missiles are currently manually allocated for naval strike 
warfare tasking. Naval Surface Warfare Center Dahlgren Division has proposed to the 
Office of Naval Research to develop automated “predesignation”’ aids that automatically 
allocate the Tomahawk Land Attack Missile, at both the Tomahawk Strike Coordinator 
level, and at the individual firing platform level. A mixed integer program is introduced 
for Tomahawk Strike Coordinator predesignation, and is implemented in alternate forms 


and tested for a variety of scenarios. 
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LIST OF TERMS AND ABBREVIATIONS 


Back-up missile 


Brilliant Anti-Tank round 
(BAT) 


Block 


Commander In Chief 
(CINC) 


Expanded Missile Identification 
(XMID) 


First Pre-Planned Waypoint 


(FPPWP) 


GREEN CROWN 


A missile, associated with a primary missile 
on another ship or submarine, that a firing 
platform is directed to align and prepare for 
launch. If the firing platform that is tasked to 
launch the primary missile experiences a 
failure of that missile, the Tomahawk Strike 
Coordinator (TSC) will order the platform 
with the back-up missile to launch. 


A weapon designed to locate and destroy 
tanks. 


A package of improvements to a missile. 
There are currently two blocks of 
Tomahawks in inventory, Block II and 
Block I. Block I improvements include 
Global Positioning System (HRS) 
navigation, lighter warheads for C variants, 


_ and enhanced inertial navigation. 


The officer tn command of a theater of 
operations. 


A series of numbers that provide detailed 
information about a missile’s software, 
including the version of operational flight 
software (OFS), and the version of HRS 
flight software (GFS). The XMID is used in 
determining position on the M’ list. 


The starting point of a Tomahawks over- 
land flight path, included in the mission 
data. Firing platforms are responsible for 
planning the over-water portion of the flight 
path: the flight from the firing platform to 
the FPPWP. 


A surface ship responsible for the air 
defense of an Amphibious Readiness Group. 
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Half-module 


Human Computer Interface 
(HCI) 


INDIGO 


Launch Sequence Plan (LSP) 


Missile Mission Matching (M?) 


Mission Distribution System 
(MDS) 


Mission ID Number (mission ID) 


A component of the Mark 41 Vertical 
Launch System (VLS) consisting of four 
individual cells. Each cell provides storage 
and launch preparation for a missile in a 
canister. Due to missile power requirements 
and VLS design, only one missile per half- 
module may be prepared for launch at a 
time. 


A tool that allows an operator to interact 
with a computer utility. 


A standard Navy message that directs a 


firing platform to align, prepare, and launch 


Tomahawks. INDIGO assigns primary, 
ready-spare, and back-up missiles. It may 
direct which cells and missiles are to be 
used, or give the number of missiles for each 
task, leaving cell and missile selection to the 


firing platform. 


A message similar to IVDIGO, but not 
granting authority to align, prepare, or 
launch missiles. It is generally issued in 
advance of an INDIGO as a planning aid for 
the firing platforms. 


A list associated with each mission ID that 
contains all of the missile. types that may be 
fired for that mission. The list is ordered in 
ascending missile capability. To maintain 
the greatest capability for follow-on tasking, 
it is desirable to select a missile from as 
close to the top of the M° list as possible. 


A tactical computer used for managing 
database information about firing platforms, 
missions and weapons, and for displaying 
TLAM mission planning information. 


A number corresponding to a single, unique, 
over-land flight path and target aimpoint. 
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National Command Authority 
(NCA) 


Naval Surface Warfare Center 
Dahlgren Division (NSWCDD) 


Primary missile 


Ready-spare missile 


RED CROWN 
Tactical Tomahawk Weapons _ 
Control System (TTWCS) 


Time on Top (TOT) 


Tomahawk Land Attack Missile 
(TLAM) 


Tomahawk Strike Coordinator 
(TSC) 


Variant 





The members of the executive branch of 
government authorized to direct the use of 
military force. 


A naval command tasked with, among other 
things, conducting research on weapon 
system development and use. 


A missile that a firing platform is directed to 
align, prepare for launch, and fire. 


A missile, associated with a primary missile 
on the same firing platform, that a firing 
platform is directed to align and prepare for 
launch. The firing platform will not fire the 
ready-spare missile unless its associated 
primary fails to launch, or fails to transition 
from the boost phase to the cruise phase of 
flight. If the primary missile fails, the ready- 
spare will be launched without further 
direction from the Tomahawk Strike 
Coordinator. 


A surface ship responsible for the air 
defense of a Carrier Battle Group. 


Proposed upgrade to the shipboard systems 
used for Tomahawk engagement planning. 


The desired time to have the tasked 
Tomahawk impact the aimpoint. 


A cruise missile, fired from surface ships 
and submarines, capable of striking land 
targets. 


The officer on the Battle Group or Naval 
Component Commander staff . 
responsible for the employment of 
Tomahawk missiles by the Battle Group. 


The type of warhead on a missile. “C” 
missiles have a unitary warhead; “D” 


AV 





missiles have multiple triple-effect 





bomblets. 
Vertical Launch System The storage, preparation, and launch system 
(VLS) for several types of weapons, including 
TLAM. 
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EXECUTIVE SUMMARY 

Employment of the Tomahawk Land Attack Missile (TLAM) in strike warfare is 
a mission area that has matured greatly since its inaugural use in Operation Desert Storm. 
Doctrine has been developed to govern the employment of TLAM in extended warfare, 
or campaign operations, as well as in measured response, or contingency, operations. 
Surprisingly, there is no official guidance governing the complex process of allocating 
missiles aboard ships and submarines to meet National Command Authority (NCA) or 
theater Commander In Chief (CINC) tasking. The planning that is crucial to the success 
of current operations, and that determines the ability of the on-scene force to conduct 
future operations, is done manually by each individual Tomahawk Strike Coordinator 
(TSC) and firing platform. No Tactical Decision Aids, Naval Warfare Publications, or 
Tactical Memoranda exist to assist the fleet operator in producing consistent, 
reproducible, and logical solutions. Forward deployed Battle Group assets will continue 
to be the on-call response force for implementing U.S. policy. Introduction of a 
systematic method to meet current tasking while maximizing capability for follow-on 
strikes, maximizing task achievement, and minimizing the impact on other operations, 
offers inestimable value to both the fleet and the nation. 

Naval Surface Warfare Dahlgren Division (NSWCDD) has proposed to the Office 
of Naval Research (ONR) to develop “predesignation” aids that will automate allocations 
for TLAM, for both the TSC, and the individual firing platform. In support of this effort, 


we examine the factors that the TSC must consider when assigning tasking, introduce an 


X1X 


allocation method for TSC predesignation, and explore alternatives for implementing the 
method. 

Through informal discussions with NSWCDD and Tactical Training Group 
Atlantic, we have determined priorities to guide the selection of firing platforms and 
missiles. In descending order of importance, these priorities are: 


a. Meet all assigned tasking. 

b. Minimize the use of firing platforms that are already occupied in other 
operations. 

c. Maximize the use of missiles from designated “expend” platforms. These are 
ships and submarines from which it is desired to use as many missiles as 
possible. This may occur, for example, if the firing platform is leaving the 
theater of operations, or is heading into an in-theater port for extensive repair 
work. 

d. Use missiles from the platforms not designated as “expend” in such a way as 
to level the number of TLAM’s remaining on these non-expend platforms. 
This is done to prevent the inventory of a particular firing platform from being 
exhausted, rendering the platform unable to perform future tasking. 

e. Ifdesired by the operator, spread tasking among as many firing platforms as 

_ possible to limit the chance of a single point of failure. 

f. Where possible, use the least capable missile for each mission. This attempts 
to limit the use of missiles with specific capabilities (e.g., GPS navigation) to 
missions that require those capabilities. 

g. Maximize follow-on firing capability. 


Using these priorities, we formulate a mixed integer linear program to perform 
allocations. This consists of a multiple-term objective finetion that expresses the relative 
priorities, subject to a series of constraints. This constrained optimization model attempts 
to allocate missiles in such a way as to reflect all of the priorities, while meeting 
operational restrictions on the use of missiles and firing platforms. 

We initially attempt to solve this mixed integer linear program as a single, 
monolithic problem. We have discovered, however, that the size of the problem renders 


this method intractable. Limits on machine precision, and the necessity to accept 








solutions that, while not optimal, are within a specified range of being optimal, result in 
solution times that are too long to be operationally useful, and are of unacceptable 
quality, i.e., the solutions returned io not adequately capture all of the priority levels. 

We next use a hierarchical restriction method in which we express the formulation 
as a series of sub-problems, one for each priority level. These sub-problems are solved in 
descending order of priority. This method accurately captures all of the priority levels, 
and returns solutions of acceptable quality; however, the solution time far exceeds what 
could be considered tactically viable. 

Our final approach is the use of a simplified heuristic, or rule-based system, to 
assist in rapidly finding solutions. The heuristic uses conditional checks to determine if 
certain sub-problems can be eliminated from consideration, uses selectively restricted 
data sets to solve some problems, and evaluates certain sub-problems using allocations | 
that were ovaned in previous sub-problems. The resulting solutions are of comparable 
quality to those obtained with the hierarchical restriction method, and drastically reduce 
the solution time required. 

We conclude that solving TSC predesignation as a monolithic mixed integer 
problem is not viable. We recommend that NSWCDD develop a complete (non-integer 
programming-based) heuristic for this application. A signal disadvantage of a heuristic is 
that it is seldom possible to objectively assess the quality of its solutions. The 
mathematical models and formal methods introduced here may be useful in 
independently assessing the quality of the heuristic solution. We further recommend that, 


at an early stage of development, this heuristic be presented to the strike warfare 





community for review. The goal of this review by decision makers is to determine the 
suitability of the heuristic for tactical systems, and to establish doctrinal ordering of 
priorities. We recommend that NSWCDD advocate the acceptance of the heuristic by the 
Naval Doctrine Command, Fleet Commanders, and Theater CINC’s, and its inclusion in 
the Naval Warfare Publications and Concepts of Operations that govern the conduct of 


strike warfare. It is crucial that strike planning be standardized in published doctrine so 


that it can be executed with the effectiveness that it deserves. 














I. INTRODUCTION 


A. BACKGROUND 


The Tomahawk Land Attack Missile (TLAM) has proven to be the weapon of 
choice for many U.S. Navy strike warfare scenarios (Figure 1). Doctrine has been 


developed to guide the employment of TLAM in extended warfare, or campaign, 


operations, as well as in measured response, or contingency, operations. 





Figure 1: An ARLEIGH BURKE Class Destroyer Launches a TLAM During 
Operation DESERT FOX [U.S. Navy, 1999a]. 


Allocating individual ships, submarines, and missiles to meet National Command 


Authority (NCA) or theater Commander In Chief (CINC) tasking is complex. Each 


individual Battle Group Tomahawk Strike Coordinator (TSC) and firing platform 








currently uses manual methods for these allocations. There are no Tactical Decision Aids, 
Naval Warfare Publications, or Tactical Memoranda to assist in producing a solution. No 
systematic way currently exists to meet current tasking while maximizing capability for 
follow-on strikes, maximizing task achievement, or representing the difference between 
the assignment of TLAM tasking to a firing shuee that is already occupied with 
another real-world mission rather than a firing platform that is more available. 

Forward deployed Battle Group assets will continue to be the on-call response 
force for implementing U.S. policy. This includes not only TLAM tasking, but diverse 
missions such as Battle Group air defense (RED CROWN), amphibious readiness group 
air defense (GREEN CROWN), non-combat evacuation of American citizens overseas, 
enforcement and monitoring of no-fly zones, humanitarian missions, maritime 
interdiction, and coalition and treaty operations and exercises. As the demands on | 
operational forces increase, fuel allocations remain constant or decrease, and weapons 
costs increase, it becomes more vital to use each ship and missile in the most vicciue 
manner. 

With the development of improved strike munitions, more options will be 
available for the NCA and CINC, such as Brilliant Anti-Tank (BAT) weapons and 
Tactical Tomahawk, allowing for rapid “call for fire” strike missions as well as the 
traditional attack against a known, stationary target. This will add complexity to weapon 


and ship allocation to meet tasking while reducing the amount of time available for strike 


planning, making manual allocation methods even less viable. 








Automated selection algorithms for the allocation of firing platforms for TLAM 
and follow-on strike munitions is desirable to provide consistent, logical, and 
reproducible solutions. Consistency is necessary because tasking may be performed at 
several levels. The TSC, or higher authority, may task to the firing platform level, leaving 
the individual firing platforms to select missiles and cells or torpedo tubes. In a more 
detailed scenario, the tasking may be specified to a particular missile. Any system of 
allocation should produce the same results, whether it is run in the Combat Information 
Center of a destroyer, or on the aircraft carrier where the Battle Group staff is embarked. 
Consistency also allows for a smooth transition between the TSC and an alternate 
commander in the event that it becomes necessary to shift TSC duties. Logical solutions 
are needed to ensure that the weapon selection follows doctrine. Reproducible solutions 
praitors effective training, and allow intelligent discrimination between alternate 
scenarios. None of these desirable characteristics exist with the current method in which 
each operator manually selects firing platforms and missiles, based on his or her 
experience and intuition. The lack of a uniform method for assigning assets to meet 
tasking denies the fleet the ability to maximize remaining power projection capability, 
ensure optimal tactical employment, and meet tasking with the least capable missile. A 
fleet-wide system of automated TLAM “predesignation,” is one way to meet all of these 


criteria. 








B. PAST WORK 
Naval Surface Warfare Center Dahlgren Division (NSWCDD) has formulated and 


implemented a mixed integer program that selects missiles to meet assigned tasking for 
an individual surface ship [Naval Surface Warfare Center Dahlgren Division, 1997]. This 
single-ship predesignation model is designed to be incorporated in the TLAM fire control 
system. The Naval Postgraduate School (NPS) thesis “Optimizing Selection of 
Tomahawk Cruise Missiles” by LT Scott Kuykendall presents an alternate formulation 
that is also capable of performing single-platform inate for a surface ship 
[Kuykendall, 1998]. Kuykendall’s mixed integer program additionally incorporates 
submarine allocation, and is extendible to a multiple firing platform Battle Group. 


C. PRESENT WORK 
Based on the success of both the NSWCDD and NPS formulations, NSWCDD 


has submitted a proposal for a multi-year project to the Office of Naval Research (ONR) 
[Fennemore, 1999c]. The goal of this project is to provide the fleet with a mathematically 
defensible solution to the problem of which specific firing platforms and missiles to use 
to meet strike tasking. Fleet users are considered both at the TSC level, and at the 
individual firing platform level. Due to differences in weapon system configurations, 
surface ships and submarines are considered separately. 

The predesignation system should consist of three separate utilities. TSC 
predesignation is designed to be embedded in the Mission Display System (MDS), and is 


intended for use by the Battle Group TSC to automate the selection of firing platforms, 


number of missiles for each platform, and, if desired, individual cells or torpedo tubes in 














order to optimize missile selection for the entire Battle Group. The output from the TSC 
solver can then be used to write LSP’s or INDIGO’s (formatted messages which provide 
advance planning information and authority to launch TLAM’s, respectively). Surface 
ship predesignation 1s designed to be included in the Tactical Tomahawk Weapon 
Control System (TTWCS), and used when the TSC does not provide tasking for specific 
cell selection. Surface ship predesignation will take an INDIGO or LSP and own-ship 
information as input, and souk a solution that optimizes half-module and cell selection 
at the individual firing platform level. Submarine predesignation is designed to be 
incorporated into a submarine’s fire control system for use in a manner similar to surface 
ship predesignation, and will automate selection of individual cells or missile and torpedo 
tube combinations. 

A brief overview of surface ship and submarine predesignation is given in this. 
thesis to provide the reader with a more complete understandin g of the entire 
predesignation project; the systems themselves are bein g developed by NSWCDD [1997, 
1999]. 

TSC predesignation will incorporate three distinct packages: the TSC solver, the 
surface ship solver, and the submarine solver. The surface ship and submarine solvers 
will be eonics of the solvers found on the individual firing platforms. The TSC solver will 
consist of three components: a Human Computer Interface (HCI), a pre-processor, and a 
solver. The HCI will allow the input and viewing of the following required data: 

a. List of targets corresponding to TLAM mission identification numbers 

(mission ID’s); 


b. Number of primary missiles to use for each mission ID; 
c. Number of ready-spare and back-up missiles to prepare for each mission ID; 





a th © 


jum 
e 


k. 


Missile Mission Matching (M?) list for each mission ID; this list specifies 
which weapon types may be used for a particular mission, and the relative 
desirability of each allowed weapon type for that mission; 

Time on Top (TOT) for each mission ID; 

Firing platform geographic location; 

Firing platform current tasking, e.g., maritime interdiction operations; 

Firing platform missile load-out data, i.e., missile block, variant, and 
Expanded Missile Identification (XMID) information down to the cell/torpedo 
room level, for each firing platform; 

Firing platform mission load-out data, i.e., list of all TLAM mission ID’s for 
which a firing platform has mission data; 

Firing platform time remaining in theater, i.e., the amount of time a particular 
ship or submarine has left in the area of operations; this is significant when an 
individual firing platform will be departing the theater earlier than the rest of 
the Battle Group; and 

The relative value of each weapon type. 


[Fennemore, 1998a] 


The operator will be able to specify tactical options, consisting of: 


a. 
b. 


Individual platforms or a group of ships to use; 

Of the platforms to be used, platforms from which it is desired to expend as 
many missiles as possible, e.g., platforms that are leaving the theater of 
operations; 

Missiles that are "reserved" for special tasking or call-for-fire missions and 
not to be considered in the solution; 

Changes to the priority order of the M? list that may be done due to direction 
from higher authority to use a missile for a task that does not require some of 
the capabilities of that missile; for example, the TSC may be directed to use 
only block III missiles for a particular task that does not require block II 
capabilities. This may occur in a contingency scenario governed by factors 
necessary to the conduct of the operation that are not of immediate tactical 
concern to the TSC; 

Geographic feasibility for each platform with respect to each launch area, i.e., 
whether or not a firing platform can attain the necessary launch position in a 
particular launch area in time to meet tasking; 

Employment penalties for each platform that assess the relative value of the 
mission that each firing platform is currently performing. For example, it may 
be more desirable to leave a cruiser on RED CROWN in the Adriatic than to 
use it for TLAM tasking in the East Mediterranean; and 





g. Launch areas in which it is desired to spread primary or back-up missiles 
across as many firing platforms as possible, or, conversely, launch areas in 
which it is desired to concentrate primary or back-up missiles on as few firing 
platforms as possible. 

[Fennemore, 1998b} 

The pre-processor will take input from the HCI and formulate a model to pass to 
the solver. The solver will be a software package capable of taking the pre-processed 
formulation and returning a solution. The solution will be processed by the HCI and 
displayed. 

The TSC solver will check tasking assigned to each firing platform down to the 
cell level. This guarantees that there is a feasible solution for all firing platforms, but does 
not address optimality from the firing platform perspective. If the TSC desires to task 
—_ to the individual missile level, the tasking proposed by the TSC solver will be 
passed to the hip and submarine solvers internal to the TSC predesignation system. 
These solvers will return half-module and cell selection solutions for each firing 
platform, which can be disseminated in INDIGO or LSP messages. If the TSC does not 
desire to task to the individual ieee level, the output from the TSC solver will not be 
passed to the internal ship and submarine solvers. Instead, the solution, aggregated to the 
firing platform level, will be sent to each firing platform as an INDIGO or LSP, where 
-» the individual surface ship and submarine solvers will be used for assignment of tasking 
to specific cells within half-modules. 


At the individual firing platform level, the optimization problem becomes less 


complex. Surface ship and submarine predesignation will consist of a single component, 








the surface ship or submarine solver, respectively. This will include an HCI, a pre- 


processor, and a solver. 


D. WEAPONS SYSTEMS 


There are currently three classes of surface ships capable of firing TLAM; USS 
TICONDEROGA (CG-47) class Aegis cruisers (with the exception of the first five ships 
in the class), USS ARLEGIH BURKE (DDG-51) class Aegis destroyers, and USS 
SPRUANCE (DD-963) class destroyers. 

All three ship types are equipped with the MK41 Vertical Launch System (VLS), 
which is the storage, preparation, and launch system for several types of weapons, 
including TLAM. The MK 41 VLS employs two launcher sizes. Full-size launchers have 
eight modules with eight cells per module, or 16 half-modules with four cells each; see 
Figure 2. Three cells of Module Five se used for an ordnance handling crane, yielding a 


full-size launcher with a total of 61 useable cells. Half-size launchers have four modules 


with eight cells per module, and also incorporate a crane, leaving 29 useable cells. 
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| Figure 2: Full-Size VLS. The full-size VLS launcher has eight modules, each with 
eight cells. The VLS ordnance handling crane occupies cells six, seven, and eight of 
module five, leaving 61 useable cells. [After U.S. Navy, 1994]. 





Maximum theoretical salvo size is determined by the total number of half- 
modules on a firing platform. This is because the MK 41 VLS only allows one missile 
per half-module to be prepared for launch at a time. In practice, the actual maximum. 
salvo size also depends on how missiles are loaded into the launcher. For a particular 
weapon type, the salvo size is the sum of half-modules that contain at least one weapon 
of that type. Figure 3 shows an example of maximum salvo calculation, using a 


hypothetical load-out. 
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Figure 3: Sample Load Plan for a Full Size Launcher. CU, CI, DU, and ~ 
refer to blocks and variants TLAM. Note that a CII missile is loaded into each 
half-module; this gives a maximum salvo size of 16 CI TLAM’s from this 
launcher. 





Cruisers (Figure 4) have two full-size VLS launchers, one forward and one aft; 
this provides a total of 122 cells and 32 half-modules. The maximum theoretical salvo 


size for a cruiser, assuming at least one TLAM in each half-module, ts 32. 
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Figure 4: USS SAN JACINTO (CG-56), a TICONDEROGA Class Cruiser [U.S. 
Navy, 1999b]. Cruisers have two full-size VLS launchers, one forward and one aft; 
this provides a total of 122 cells and 32 half-modules. The maximum theoretical 
salvo size for a cruiser, assuming at least one TLAM in each half-module, is 32. 


ARLEIGH BURKE class destroyers (Figure 5) have a half-size VLS launcher 


forward, and a full-size launcher aft, providing a total of 24 half-modules and 90 cells, 


* 


and a maximum salvo of 24. 
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Figure 5: USS MITSCHER (DDG-57), an ARLEIGH BURKE Destroyer [U.S. 
Navy, 1999c]. ARLEIGH BURKE class destroyers have a half size VLS launcher 
forward, and a full size launcher aft, providing a total of 24 half-modules and 90 
cells, and a maximum salvo of 24. 

SPRUANCE class destroyers (Figure 6) have a single full-size launcher forward, 
16 half-modules and 61 cells, for a maximum salvo of 16. Additionally, due to ship 
combat systems design, SPRUANCE class destroyers only carry TLAM’s or Vertical 
Launch Anti-Submarine Rockets, a rocket thrown torpedo, in their VLS launchers, while 
cruisers and ARLEIGH BURKE class destroyers can additionally load SM-2 variant 


surface-to-air missiles. Therefore, in general, SPRUANCE class destroyers will be 


loaded out with a higher proportion of TLAM’s than ships of the other two classes. 
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— Il. MODELING APPROACHES 


A. PROBLEM DEFINITION 

The scope of this thesis is the formulation of the mixed integer program that will 
be the template for the TSC solver, and the implementation of the formulation in GAMS 
[GAMS, 1999]. Per discussions with NSWCDD, the formulation accounts for multiple 
Battle Groups, launch areas and time periods. The objective function minimizes penalties 
for failure to meet tasking, employment penalties, and penalties for not using the least 
capable missile, while maximizing follow-on salvo capability and allocating missiles 
across firing platforms in accordance with the tactical decisions of the TSC. 


B. ASSUMPTIONS 


The TSC Predesignation model is a mixed integer program with a linear objective 
function and linear constraints. Several assumptions and modeling approaches in the 
formulation need explanations. 

Considerations that are unique to submarines are not addressed. NSWCDD has 
determined that it is sufficient to demonstrate the TSC model using only surface ships as 
firing platforms [Fennemore, 1999a]. 

An index is not explicitly defined to differentiate between firing platforms of 
different Battle Groups, as any Battle Group asset can be uniquely identified by its hull 
number. For example, USS MITSCHER has the hull number DDG-S7. If USS 
MITSCHER is assigned to the ENTERPRISE Battle Group, referring to DDG-57 


conveys the firing platform as well as the associated Battle Group. 
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The model treats tasking as a series of tasks and parts. Each task corresponds to a 
single TLAM mission ID, and consists of three parts. Parts 1, 2, and 3 correspond to the 
primary, ready-spare, and back-up missile, respectively. While each task corresponds to 
one mission ID, a mission ID may have multiple associated tasks. Each part of each task 
can require the firing of at most one missile. To calculate the number of missiles 
remaining after a strike, it is assumed that primary missiles will be launched, but ready- 
spare and back-up missiles will not. 

Tasks are associated with launch areas. A launch area is a geographic location 
containing any positive number of First Pre-Planned Way Points (FPPWP’s). A firing 
platform in the launch area can fire missiles for any mission with a FPPWP in that launch 
area, and cannot fire missiles for any mission outside of that launch area. This is due to 
the maximum ranges of the various missile types, and corresponding restrictions on the 
length of the over-water portion of the flight path. For example, a ship in the Adriatic 
launch area can fire missiles for any mission with a FPPWP in that launch area, but 
cannot fire any missile for missions that have FPPWP’s in the North Arabian Gulf. 

Missile spreading, i.e., how missiles are chosen from several operationally 
feasible firing platforms, depends on three features. First, it may be desirable to use as 
many missiles from a particular platform as possible. This situation may occur, for 
example, because a ship is leaving the theater of operations, or is heading into an in- 
theater port for extensive repair work. Such a platform will be designated in the 
formulation as an “expend” platform, and, with all other factors being equal, the mode] 


will assign to it the maximum possible tasking. The next priority in missile spreading is 
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to level assets across the non-expend platforms. The model will attempt to leave all non- 
expend firing platforms with the same number of missiles; this prevents depleting the 
magazine of any non-expend firing platform. Finally, as a tactical concern, it may be 
desired to spread primary and/or back-up missiles across as many firing platforms within 
a launch area as possible, or, conversely, to concentrate assigned tasking on as few 
platforms as possible. This accounts for the trade-off between minimizing the potential 
damage of a single point failure, and minimizing the number of platforms that are 
required for TLAM tasking. The option to spread task parts individually (i.e., primary 
missiles but not back-ups, or vice versa) is allowed; this gives the TSC flexibility in cases 
where a primary to back-up ratio other than 1:1 1s desired. Such a case may arise, for 
example, if a number of back-ups are tasked that have no associated primary missile. This 
form of tasking is used to enable the rapid execution of strikes that must be planned for, . 


but for which the TSC has not yet received launch authority. 


C. MODEL PRESENTATION 


The data in the formulation is derived from the data and tactical options 
enumerated in the description of the overall TSC predesignation system. The primary 
decision variables for the formulation are binary variables, indexed over firing platforms, 
half-modules, cells, tasks, and task parts (primary, ready-spare, or back-up). These 
variables represent the selection, or non-selection, of a missile on board firing platform /, 
in half-module h/ and cell c, to meet task t, part p. 

The objective function is presented in order of solution priorities. Priorities are 


assigned as follows, from highest to lowest: 
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a. Meet all assigned tasking through the use of a penalty assessed for any unmet 
task parts. Primary missiles are considered the most important to select, 
followed by back-up missiles, and then ready-spares; 

b. Minimize employment penalties, which quantify the opportunity cost of using 
a platform for TLAM operations instead of other tasking; 

c. Maximize missiles allocated from designated expend platforms; 

d. Level the number of missiles remaining on non-expend platforms; 

e. Spread primary missiles among as many firing platforms as possible, in 
launch areas so designated; 

f. Spread back-up missiles among as many firing platforms as possible, in 
launch areas so designated; 


g. Minimize the M? list position of all selected missiles, ensuring that, all other 
priorities being satisfied, the least capable missile is allocated for each 
mission; and, finally, 

h. Maximize residual salvo capability. 

Although the ordering of the priorities presented in this thesis was developed 
jointly with NSWCDD [Fennemore, 1998c], and was informally reviewed by members of 
Tactical Training Group Atlantic [Williams, 1998], it does not reflect any official 
guidance. 

The constraints assign values to the accounting variables used in the objective 
function and ensure that a solution is operationally feasible. General constraints are 
imposed to govern the way in which missiles may be chosen for tasking, and ships used 
to handle tasking requirements. Noteworthy operational requirements are as follows: 

_ (a) Only one cell per half-module is allowed to be selected; this is due to a design 
characteristic of the Vertical Launch System (VLS); (b) Ready-spare missiles for a task 
can only be assigned to the firing platform to which the primary missile for that task has 
been assigned; and (c) back-up missiles can only be assigned to a firing platform other 


than the one to which the primary missile for that task has been assigned. The general 


format of several constraint sets is derived from Kuykendall [1998]. 
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A more specific description of each term of the objective function and each 


constraint follows the formulation. 
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FORMULATION 
Indices: 

Ww weapon type loaded in cell, including block and variant (e.g., CII, DI). 

f firing platform (e.g., DDG-57, CG-73, DD-997). 

h half-module, number is dependent on type of ship (e.g., hl-h24 for DDG- 
57). | 

Cc cell, each half-module contains four cells (c1-c4) or (c5-c8). 

t task number (e.g., tl, t2, ...); each task number corresponds to a mission 
ID number. Mission ID numbers may have multiple associated tasks. 

P task part, each task may consist of a primary part, a ready-spare part, and a 
back-up part indexed as follows: 1 = primary, 2 = ready-spare, 3 = back- 
up. 

a launch area, i.e., geographic location with multiple FPPWP in close 
proximity (e.g., N. RED SEA, EAST MED). 

i instance of a set of conflicting tasks (e.g., {t1, t2, t3}) 

Set: 

CON; set of tasks that are to be executed within a 
period of time such that all tasks in the set - 
conflict, i.e., only one task per set can be 
assigned to any one half-module, where con; 
denotes a set of conflicting tasks. 

Binary Data: 

inArear, equals 1 if task t is in launch area a, 0 
otherwise 

geoFes, equals 1 if platform f can attain firing 
position for task t, 0 otherwise 

numRedtp equals 1 if task t, part p requires a missile, 0 
otherwise 

expends ~ equals 1 if it is desired to expend as many 
missiles as possible from platform f, 0 
otherwise 
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priSpread, 
buSpread, 


loadwmc 


empPens 


Weighting Data: 


mCubePrip 


Numerical Data: 


levelPlat 


mCubePos,; 


value, 





equals 1 if it is desired to spread primary 
missiles across as many platforms in launch 
area a as possible, 0 otherwise 


equals 1 if it is desired to spread back-up 
missiles across as many platforms in launch 
area a as possible, 0 otherwise 


equals 1 if a weapon of type w is loaded in 
location (f,h,c), 0 otherwise 


equals 1 if it is desired to have firing 
platform f continue its current mission 


' instead of being used for TLAM tasking, 0 
otherwise 


- weight of M? list position for task part p 


total number of firing platforms in the 
operational theater minus the number of 
“expend platforms”, 1.e., the group of 
platforms across which it is desired to level 
the number of missiles after a strike 


M? list position of weapon type w for task ¢ 
if the M° list for task ¢ contains weapon type 
w, O otherwise 


relative value of weapon type w (e.g., CII 
are more valuable than CI, which are more 
valuable than DID) 
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Derived Data: 


availfip 
geoFe 
mCubePos; 
loady 
inArea, 


inArea, 


Binary Variables: 


Xfhetp 
UNABLE, 
SHIPREQnf 


NUMPLATPRg 


NUMPLATBUs 





equals 1 if firing platform f, half-module h, 
cell c contains a TLAM available for 
assignment to meet part p of task f, 0 
otherwise 

{(h,c) | availmctp = 1} 

{(f,t) | geoFeg = 1} 

{w | mCubePos,; =1} 

{(f,h,c) | loadyee = 1} 

{a linArea,, = 1} 


{t linArea,, = 1} 


equals 1 if weapon in location (f,h,c) is 
selected for task t, part p, 0 otherwise 


equals 1 if missile cannot be allocated for 
task t, part p, 0 otherwise. 


equals 1 if platform fis required to be in | 
area a for firing, 0 otherwise 


equals 1 if one or more weapons of type w 
remain on platform fin half-module h after 
firing all primary missiles, 0 otherwise 


equals 1 if platform fis assigned one or 
more primary task parts, 0 otherwise 


equals 1 if platform fis assigned one or 
more back-up task parts, 0 otherwise 
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Counting Variables: 


NUMEXPEND 


NUMREMAIN¢ 


AVGREMAIN 


DIFFMEAN;s 


MCUBESUM, 





number of missiles selected from platforms 
designated to receive as much tasking as 
possible 


number of missiles on platform that are 
not selected for a primary task part from 
non-expend platforms 


mean number of missiles of all types on 
board “non-expend” firing platforms in 
theater, after all missiles have been selected 
for the given strike 


difference between the residual number of 
missiles on platform f and AVGREMAIN 


sum of M? list positions for missiles that are 
selected for part p 
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Formulation: 


MINIMIZE 
la) } UNABLE, 


ip 


1b) + >) (empPen, * SHIPREQ,, ) 
fa 

ic) - NUMEXPEND 

1d) + }’ DIFFMEAN, 
f 


le) - ). NUMPLATPR, 
f 


if) - ), NUMPLATBU, 
f 


1g) + } mCubePri, * MCUBESUM, 
: | 


1h) - } (value, * SALVO,,,) 
wth 


Subject to: 


2) » > x < numberReq, _“V (t,p)|numReq,, =1 


fhctp —™ 
f lgeoFe, =1 helavailp.., =! 


3) Xpurrr SSHIPREQ, V(f,t)lgeoFe=1 
a linArea, =1 
(h,c)lavail,.,. =1 
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4)X que <SHIPREQ,  V(f,t)l geoFe =1 
a linArea, =1 


(h,c) lavail,.,.=1 


5) );SHIPREQ,, <1 Vf 


() » > Xaeers = NUMEXPEND 


f lexpend, =1 hetlavaila.)- =1 


7) Sloadyme- > Xmery = NUMREMAIN, Vflexpend, =0 


whe het lavaile.y- =1 


a  NUMREMAIN, = AVGREMAIN 
levelPlat 4 


9) (NUMREMAIN, - AVGREMAIN)*(1-expend,) 2 -DIFFMEAN, Vf lexpend, =0 
10) (NUMREMAIN, - AVGREMAIN) * (1-expend,) < DIFFMEAN, Vf lexpend, =0 


1) YY YY Xpeqe 2 NUMPLATPR, VE 


al priSpread, =1 tlinArea, =1 hclavailp.4-}- =! 
12) NUMPLATPR, <1 Vf 


3) > » SY Xecesx 2 NUMPLATBU, Vf 


a! buSpread, =1 tlinArea, =] hclavaily.-3- =1 


14) NUMPLATBU, <1 Vf 


5 » > . Xm = MCUBESUM, Vp 


wimCubePos, =! fhcfload, =1 tlavailp..) =! 


16) > loadsa, - > Xaep 2SALVO,, V(w.f,b)l >} load,,, 21 


fhetp “~ 
Clavailaaey =! ctlavailaary: =1 
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17) SALVO,, $1 V(f,h)l } load, 21 


cw 


18) SM Xaww 2 Xa  V Et) geoFe, =1 


-helavail gery. =1 he! avail pegs =1 


19) SY Xam FY Xp $1 V Eb |DICY Toad, 21 VU 


clavailg ae =1 Clavailp eg =1 cw 


t | (numReq,.,. +numReq ,») = 2) 


20) YY Xaw + SM Xa $1 V E,t)lgeoFe, =1 


he lavatl pyermy = | helavail parse =] 


21) Xa, $1 V(h,c)l}load,,, 21 


tplavailaa, =! 


22) Xm $1 V (tp)!numReq,, =1 


fhe! avail grep =] 


23) SY YS Kau $1 VEb)I load. = Vi 


cp lavaila, =1 te con; 
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NOTES 
Objective function explanation: 
la) Expresses the penalty for failure to meet tasking requirements, 


lb) plus the sum of employment penalties of platforms being used for Tomahawk 
tasking instead of continuing the mission they were performing, 


Ic) minus the benefit for the use of missiles from designated “expend” platforms, 
ld) plus a penalty for differences between the number of missiles remaining on non- 
expend platforms after tasking, and the average number of missiles remaining on 


non-expend platforms after tasking, 


le) minus the benefit of spreading primary missiles across as many platforms as 
possible, where it is desired to do so, 


-1f) minus the benefit of spreading back-up missiles across as many platforms as 
possible, where it is desired to do so, 


lg) plus a penalty that increases with the M? list position of the selected missile, and 
[h) minus the total residual salvo capability for each weapon type. 
Constraint explanations: 


2) The sum of missiles selected and the number of missiles that cannot be selected is 


equal to the number of missiles required for each task and part. The variable Xfrep is 


only “counted” for missile selection if, for a specified launch area, the task is in the 
area, and the firing platform can attain launch position in the area. 


3) A firing platform f is required to be in area a if a missile on that firing platform is 
selected for a primary task in area a. 


4) A firing platform / is required to be in area a if a missile on that firing platform is 
selected for a back-up task in area a. 


5) A platform f can launch from no more than one area. 


6) The variable NUMEXPEND; equals the sum of all missiles selected from expend 
platforms. | 
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7) The variable NUMREMAIN; equals the initial load-out of missiles minus the sum of 
all primary missiles fired, for each platform that is not designated expend. This 
assumes that only primary missiles are fired, and that all primary missiles are fired. 


8) The average number of missiles remaining on non-expend platforms across all 
platforms in theater is equal to the sum of NUMREMAIN divided by the total number 
of non-expend platforms in theater. 


9-10) The variable DIFFMEAN; is set for each “nonexpend” firing platform by 
comparing the residual number of missiles to the average number of remaining 


missiles. 


11-12) The variable NUMPLATPR; equals one for each firing platform on which 
primary missiles are selected, for launch areas that require primary spreading. 


13-14) The variable NUMPLATBU; equals one for each firing platform on which 
back-up missiles are selected, for launch areas that require back-up spreading. 


15) The variable MCUBESUM equals the sum of M’ list positions for all missiles that 
are selected. 


16) The variable SALVOvym is restricted to equal zero if all missiles in a half-module 
have been expended, for each firing platform and weapon type. 


17) The salvo size for a weapon type for each firing platform and half-module is either 
0 or 1. This is due to half-module power-up constraints — a surface ship may only 
apply power to one cell per half-module. The formulation will only credit the 
follow-on salvo capability with at most one weapon of each type per half-module. 


18) Ready-spare missiles for a task must be assigned to the same firing platform as the 
primary missile for that task, for each firing platform and task. 


19) Ready-spare missiles for a task cannot be assigned to the same half-module as the 
primary missile for that task, for each firing platform and task. 


20) Back-up missiles for a task must be assigned to a firing platform other than that 
assigned the primary missile for that task, for each firing platform and task. 


21) Only one missile may be selected per half-module on a given firing platform. 


22) Only one missile may be selected per task part. 
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23) Only one missile per set of conflicting tasks may be selected per half-module for 
each firing platform. 
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D. IMPLEMENTATION 


The TSC predesignation model is implemented in GAMS version 2.50 [GAMS, 
1999] using the CPLEX solver version 6.5 [ILOG, 1999], and run 0 a 333 MHz PC with 
a Pentium II processor. Formulation alternatives were tested to improve solution time and 
the quality of the solution (e.g., within 5% - 10% of a provably optimal solution). 

l. Weighted Objective Function Implementation 

Initially, the model was implemented as a monolithic problem, with a single 
objective function containing all terms. Each term (e.g., using missiles from expend 
ciaeeaenna or maximizing residual salvo capability) represents a distinct priority to be 
considered in allocating weapons and firing platforms. Numerical multipliers, or weights, 
preserve the ordering of the priorities of the solution. For example, the highest priority of 
the formulation is to meet all assigned tasking. The weight applied to this objective 
function term is of sufficient magnitude to ensure that no combination of allocations 
achieved with lesser priorities can outweigh an allocation with the least amount of ar 
tasking possible. 

The most serious problem with using objective function weights to enforce 
prioritized solutions is his of inaccuracy resulting from non-zero integrality gaps (Le., 
the ratio between the best integer solution found when the algorithm terminates and the 
bound on the best possible solution). An example of the process of determining weights 
is presented to illustrate how this inaccuracy occurs. 

In the TSC predesignation formulation, the lowest priority is residual salvo 


capability. Each missile of a specific type is assigned a relative value. Because this is the 
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lowest priority, the numbers can be small, such as assigning a CI missile a value of one, 
and a CII missile a value of two. This allows for a sufficient number of less capable 
missiles to be considered more valuable than a small number of more capable missiles; 
this is an intentional effect of the formulation, and places a premium on overall striking 
power rather than the ability to accomplish a specific set of missions. With some 
assumptions about the number of firing platforms under consideration, and the expected 
number of TLAM’s available, the maximum possible value for the SALVO variable can 
be determined. Using the assumptions specific to the scenarios run in this thesis, 
discussed in detail in Chapter IJ, the maximum possible value of SALVO is: 


2, (value, \(number, ) 


where w is the type of missile, value,, is ite value of a type w 
missile, and number, is the number of missiles of type w 


= (2)(104) + (1)(78) 
=280. 

The next priority is the sum of all M? list positions for the selected missiles. This 
priority is further separated into task parts, so that primary missile M list positions are 
considered more important than back-ups, and back-ups are more important than ready- . 
spares. Using the maximum possible value for residual salvo capability of 280, a weight 
which is slightly larger, 300, is applied to each M? list position of missiles selected for 
ready-spare missions. This weight ensures that residual salvo considerations will not 
affect the selection of missiles to minimize M? list position. Assumptions about the 


maximum possible tasking, and the M? list for tasks gives a worst case scenario of 64 
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ready-spare missiles being selected, all with M? list position two. The maximum possible 


total of the sum of M° positions for ready-spares is then: 


where M° Weight is the weight for M? list positions of missiles selected 
for ready-spares, j is an index of M? list positions (e.g., one or two), and 


M? Weight » (j* number, ) 

number; is the number of missiles selected that have M? list position j 

= (300)(2)(64) 

= 38,400. 
This number then influences the selection of the weight on the objective function term for 
M? list positions for back-up missiles. This weight must account for the sum of all 
previous worst case scenarios, and so must be larger than 38,400 + 280, or 38,680. The 
remaining weights are calculated in a similar manner. For the scenarios used in the thesis, 
this results in the highest priority weight, i.e., the penalty for unmet primary tasking, of 
8.32x10"” | 

This becomes problematic when the solver is trying to differentiate between 

solutions. In order to allow the solver to find a solution in an amount of time which 
represents the desired balance between the required quality of the solution and the time 
needed to return an optimal solution within the limits of machine precision, a relative 
integrality gap is set; this gap permits the solver to return a solution which, while not 
optimal, is within the specified tolerance of being optimal. A mode] that has a wide range 
of objective function weights is in danger of having the lower weights largely ignored by 
the solver; with a relative gap greater than 1.0x107°%, an objective function value of 


8.32x10"? is indistinguishable from a value of 8.32x10'° + 280. Using the formulation 
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presented in this thesis, a small scenario assigning 29 missiles to four firing platforms 
with a combined salvo capability of 104 TLAM’s, an allowed relative integrality gap of 
5% led to a solution with an spesite integrality gap of 5.5x10'°. The effect of this gap 
was seen in the variables for the lower priority objective function terms. The model 
selected missiles that resulted in a sub-optimal solution for considerations such as the 
sum of M’ positions and residual salvo capability. 

An additional difficulty with the method of weighted objective functions is one of 
numerical stability. Every computer has limitations on the precision it maintains in 
memory while performing arithmetic operations. The computation necessary to execute 
the branch and bound algorithm on problem instances with numbers (i.e., weights) 
differentiated by large orders of magnitude arguably approaches the limits of precision of 
the computer used. 


De Hierarchical Restriction Implementation 


Based on the poor performance of the weighting, we have abandoned weighting 
in favor of multiple objective function hierarchical restriction. The monolithic objective 
function is separated into eight sub-problems, each containing a single-term objective 
function; these sub-problems are solved in decreasing order of objective function priority. 
Optimal values for the objective function term variables obtained in one sub-problem are 
fixed, and provided to the subsequent sub-problem as input. By allowing the primary 
variable, individual missile selection, to change while maintaining the values of objective 
function term variables set in previous sub-problems, it is possible to solve for each 


priority, ensuring that no lesser consideration can impact a more important concern. For 
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example, the highest priority of the model is to meet all assigned tasking. Using the 
hierarchical restriction formulation, we solve the first sub-problem to minimize failure to 
meet tasking for each task and part. Once the optimal solution is found, the values for 
these variables are set. These values are passed to the next sub-problem as data. This 
model, which minimizes any employment penalties, is then solved, with a constraint that 
unmet tasking is fixed according to the solution from the first sub-problem. Employment 
penalties are then fixed to the values returned by the second sub-problem. The next sub- 
problem, which maximizes the allocation of missiles to any expend platforms, receives 
the unmet tasking and employment penalty data, and is constrained to maintain these 
values. The process is continued until all sub-problems have been optimized. 
Hierarchical restriction also simplifies the problem by allowing us to consider a 
restricted subset of the constraints. Specifically, with a single objective function term, the 
only constraints that need to be considered are those that relate to the objective function 
variables under consideration, those that fix variables whose values have been set in a 
previous sub-problem, and constraints that enforce edtrictions on missile selection. This 
provides a significant reduction in the number of rows and columns generated for each 
model. The hierarchical restriction sequence is summarized in Table 1; the objective 


function term and constraint numbers are consistent with those given in the formulation. 
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HIERARCHICAL RESTRICTION SEQUENCE 


Sub-problem: Variable to be fixed: 


1 Minimize la 2, 3,4, 5, 18, 19, 20, 21,22, | UNABLE,, 
23 
23 
23 

4 Minimize 1d 2,3,4, 6,7, 8,9, 10, 18,19, | DIFFMEAN; 
20, 21, 22, 23 

5 Maximize le 2, 3, 4, 6, 7, 8,9, 10,11,12, | NUMPLATPR; 

| 18, 19, 20, 21, 22, 23 
Maximize If 2, 3, 4, 6, 7, 8, 9, 10, 11, 13, NUMPLATBU; 

14, 18, 19, 20, 21, 22, 23 


15, 18, 19, 20, 21, 22, 23 

Maximize ih 2, 3,4, 6, 7, 8,9, 10, 11,13, | SALVOve, 
15, 16, 17, 18, 19, 20, 21, 22, 
23 


Table 1: Hierarchical Restriction Sequence. Sub-problems are numbered in 
decreasing order of priority. The objective function and constraint numbers refer 
to the formulation. Each sub-problem is solved and the variables that appear in 
the objective function are fixed. Constraints are added to each sub-problem to 
ensure that allocations are operationally feasible. 





Hierarchical restriction allows us to more accurately capture priorities at all levels 
given hardware and software limitations. However, our mmplementation of the method- 
does not necessarily provide optimal solutions. By fixing the values of the relevant 
variables after the solution of each sub-problem, we are precluding solutions in 
subsequent sub-problems which may result from an alternate optimal soliton in a 
previously considered sub-problem. This difficulty can be avoided by constraining the 
objective function value derived from solving a given sub-problem to be (almost) as good 
as Its value when solving a subsequent sub-problem (e.g., [Steuer, 1986], pp. 292-296). 
For example, if solving sub-problem 1 yields three primary task parts that cannot be met, 


with our method, the value of the unmet tasking variable is set equal to one for the three 
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task parts returned in the optimal solution. Alternately, we could add a constraint 
ensuring that among all tasks, three primary missiles cannot be allocated. The drawback 
to the latter approach manifests itself in resulting sub-problems with a greater number of 
variables and constraints than with the sub-optimal method we used. In fact, the size of 
our problem rendered the optimal method unusable. Even with our implementation, as the 
size of the scenario increases, there is a marked increase in solution time. This is 
presented in Chapter II, but the solution times for even moderately sized scenarios are 
poor enough to mandate further implementation refinements. We next present a heuristic 
variation of our hierarchical restriction procedure to improve solution times, albeit at the 
expense of some solution quality. 


3. Heuristic Structure 


The hierarchical restriction implementation results in exceedingly long solution 
times for certain sub-problems in various problem scenarios (e.g., on the order of 2,000 
seconds for a single sub-problem in some cases). Based on observations of the behavior 
of the model, we implement a heuristic procedure in conjunction with the hierarchical 
restriction methodology. The logic for the heuristic is based on four ideas: 

First, for several of the sub-problems, it may be possible to decrease solution time 
by solving a sub-problem multiple times, each time using a restricted set of the data. 
Specifically, solving successively for individual task parts (i.e., primary, ready-spare, and 
back-up) rather than solving for all task parts simultaneously reduces solution times for 
some of the sub-problems by reducing the sub-problem size. Solving for separate task 


parts is used for sub-problems one and seven, which address unmet tasking and the sum 
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of M? list positions respectively. For sub-problem one, if all required allocations can be 
made using this method, then the unmet tasking variables are fixed to zero, the allocation 
variables are released, and the model continues to the next sub-problem. For sub-problem 
ele if a feasible solution is found for all tasking with the sequentially fixed allocation, 
then the M? list position variables are fixed, the allocations released, and the model 
proceeds to the next sub-problem. In both cases, if the required tasking cannot be met 
with this method, the sub-problem is re-solved for all task parts together, 1.e., with the 
full-scale hierarchical restriction implementation. 

Second, obtaining a solution using only a subset of data in one sub-problem may 
enable us to omit a subsequent sub-problem entirely. This method is used in instances of 
sub-problem three in which there are firing platforms that are ascioned an employment 
penalty. If all tasking can be allocated in the first sub-problem without using any of these 
firing platforms, then the value for employment penalties 1s set to zero, and the sub- 
problem for employment penalties is skipped. 

Third, we observe that for some sub-problems, almost all of the solution time 1s 
spent trying to find a feasible solution. If this initial feasible solution is not within the 
specified optimality tolerance, the solver finds one that meets the tolerance within a 
matter of seconds. Relatively few constraints are added to each successive sub-problem. 
Additionally, all operational constraints are met in each sub-problem (thus ensuring that 
any feasible allocation determined from a previous sub-problem is also feasible for a 
subsequent sub-problem). Thus, the allocation from the previous sub-problem may 


provide a relatively good solution to the subsequent sub-problem. 
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This suggests the following method for reducing solution time. Rather than 
solving a sub-problem requiring a long solution time and whose objective function value 
may not differ substantially if it is evaluated with the allocation variables fixed from the 
previous sub-problem, we proceed as follows: Fix the allocation variables to the values 
from the preceding sub-problem. Next, compute the values for the corresponding 
objective function variables. The values for allocation variables are then released, and the 
next sub-problem is solved. This results in a significant decrease in computation time, 
because the solver is performing a simple arithmetic computation rather than a potentially 
large mixed integer optimization. 

The method of fixing allocations after a preceding sub-problem and computing 
the objective function value rather than optimizing is used for sub-problems five and six, 
corresponding to the spreading of primary and back-up missiles, and sub-problem eight, 
which maximizes residual salvo capability. 

Finally, a simple if-then check may be used to determine the need to solve sub- 
problems. If a sub-problem is associated with an operator selection, and the operator does 
not desire to perform the action in question, then the sub-problem does not need to be 
run. This method is employed with sub-problem three, which maximizes allocations to 
expend platforms. If no expend platforms are designated, the sub-problem is omitted. 
This method could also be applied to sub-problems five and six; however, the time saved 
by omitting the sub-problem rather than computing the objective function value for a 


fixed allocation is negligible. 
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The structure of the heuristic implementation is presented in Figure 7. The 
conditional and procedural statements reflect the optimization of a particular sub- 
problem, whether the missile allocation is treated as fixed or variable, and the necessity 
of solving a particular sub-problem. The task blocks for solving sub-problems one and 
seven are simplified; they do not represent the process of solving for each task part 
separately. pine: 8 and 9 give detailed procedures of the heuristic method for sub- 
problems one and seven, respectively. 

There are several caveats associated with the use of this heuristic. Efficiency is 
highly problem-specific. In particular, the use of restricted data sets does not ensure a 
faster solution time; in fact, it is theoretically possible . construct cases for which the 
solution time using the heuristic is longer than that of the hierarchical restriction 
implementation, if the restricted data sets do not render a feasible solution. For example, 
the solver may spend time determining that tasking eas be met without using a 
platform with an employment penalty, and the entire sub-problem must be re-run 
considering the additional firing platforms. 

Additionally, the solution obtained using the heuristic will be no better than the 
hierarchical restriction solution, and, in some cases, may be worse. This is because 
allocations are fixed in a greedy and myopic manner. For example, in solving for M° list 
position by task part, it is possible that primary allocations will be made that will 
adversely impact the choice of ready-spare missiles with respect to M? list position, and 
the combination of fixed primary and ready-spare allocations will similarly adversely 


affect the selection of back-up missiles. Further, the decision to temporarily fix all 
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Figure 7: Heuristic Solution Method. Notation is consistent with the formulation. 
See Figures 8 and 9 for expanded views of sub-problems I1 and 7. 
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Figure 8: Expanded Heuristic for Sub-Problem 1 to Minimize Unmet Tasking. 
See Figure 7 for complete heuristic flow. | 
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Solve sub-problem 7 
for primary task parts. 
Fix MCUBESUM++». 
Fix primary allocation. 
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for ready-spare task 
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Figure 9: Expanded Heuristic for Sub-Problem 7 to Minimize M? List Position. 
See Figure 7 for complete heuristic flow. 
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allocations before calculating the number of platforms primary and back-up missiles are 
spread across, and residual salvo capability is one that implicitly accepts a sub-optimal 
solution in favor of greatly reduced solution times. 

Results for scenarios run using the heuristic are given in Chapter I. Solution 
quality appears to be acceptable for all cases tested. Note, however, that the heuristic has 
been designed expressly for these particular scenarios, and additional procedures (e.g., 
subroutines and condition cece) would be required before it could be applied to 


scenarios differing substantially from those presented here. 
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II. RESULTS 


A. SCENARIO DESCRIPTIONS 


Scenarios for demonstrating the TSC preieaanalion model are based on data 
provided by NSWCDD [Fennemore, 1999b]. There are six “base” scenarios, which are 
designated numerically (i.e., scenarios 1 — 6). The base scenarios show the effectiveness 
of the model as the amount of tasking increases. Five additional scenarios, denoted alpha- 
numerically (e.g., 1A), are presented to exercise all the options that the formulation 
allows. The additional, or modified, scenarios have the same amount of taskin g as the 
associated base scenario, but demonstrate allocations resulting from different operator 
input. For example, Scenario 1 requires 13 primary missiles, 8 ready-spare missiles, and 
8 back-up missiles. Scenario 1A requires the same amount of tasking, but also designates 
an expend platform, which is not required in Scenario 1. 

Scenario 6 introduces the concept of “fixed tasking”. This assumes that the 
tasking from Scenario 5 has already been assigned to firing platforms and missiles have 
been selected to meet this tasking. Additional tasking is then allocated, with the 
restriction that the previous allocation must remain unchanged. This capability is 
demonstrated at the request of NSWCDD, and represents a first step towards a robust 
predesignation system that is capable of performing allocations under varying operational 
conditions. 


The scenarios are summarized in Table 2. 
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Table 2: on So Ty 


Notes: 1. CG is used to denote a TICONDEROGA class cruiser. DDG indicates an 


ARLEIGH BURKE class destroyer, and DD refers to a SPRUANCE class 
destroyer. 

2. PRI, RS, and BU represent primary, ready-spare, and back-up task parts, 
respectively. 

3. Data sets expends and empPen, refer to the formulation. 

4. Scenario 1A demonstrates the use of an expend platform. 

5. Scenario 1B illustrates the spreading of back-up missiles among as many firing 
platforms as possible. 

6. Scenario 2A has two launch areas. One CG and one DD can only perform 
tasking in one of these areas; one CG and one DDG may only be assigned tasking 
in the other area. 

7. Scenarios 5A and 5B show the effects of expend platforms and employment 
penalties, respectively, with larger sets of tasking. 

8. In Scenarios 5, 5A, 5B, and 6, thirteen of the tasks may only be allocated to 
three specific firing platforms. This illustrates the possibility of firing platforms 
not having mission data or the correct missile type for certain missions. 
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For each scenario, tasks are given in groups of thirteen. Each task group requires 
thirteen primary, eight ready-spare, and eight back-up missiles. Using an assumed 
preparation time for a TLAM of approximately 45 minutes, the launch times for tasks are 
assigned so that all thirteen tasks in a group “conflict” with every other task in that group. 
This means that, because the missiles are being prepared for launch at the same time, 
only one missile per conflict set may be assigned to a half-module. Additionally, the last 
three tasks in each group conflict with the first ten tasks of the following group. For 
example, Scenario 2 assigns 26 tasks, in two groups of thirteen. Task numbers one 
through ten have a notional launch time of 1100 on the day of the scenario; task numbers 
eleven through thirteen have a launch time of 1130. Tasks fourteen through twenty-three 
are to be launched at 1200, and the final three tasks each have a launch time of 1230. 
With the assumed preparation time of 45 minutes prior to launch, this scenario generates 
the following conflict sets: Tasks 1-13, Tasks 11-23, and Tasks 14-26. Figure 10 gives a 
graphical representation of how conflict sets are formed. The span of time for each set of 
tasks corresponds to the 45 minutes of preparation time. Each preparation period ends at 


the assigned launch time for that particular group of tasks. 
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(---Tasks 1-10 
(----Tasks 11-13 
(---Tasks 14-23 


(---Tasks 24-26---) 


Time 1000 1030 1100 1130 1200 1230 


Conflict Sets: C1 = {Tasks 1-10, 11-13} 
C2 = {Tasks 11-13, 14-23} 
C3 = {Tasks 14-23, 24-26} 





Figure 10: Conflict Sets. Tasks 1 through 10 are all assigned a launch time of 1100. 
Tasks 11 through 13 are to be launched at 1130, tasks 14 through 23 at 1200, and 
tasks 24 through 26 at 1230. Only one missile at a time may be prepared for launch 
in any one half-module. Using an assumed 45 minute preparation time for each 
missile prior to launch, this results in three sets of “conflicting” tasks; only one 
missile per conflict set may be assigned to any half-module. [After Naval Surface 
Warfare Center Dahlgren Division, 1997] 

Three of the eight task groups require CII missiles; all other task groups allow 
either CII or CII missiles to be used. For the tasks requiring CII missiles, the M list 
contains.a single entry, assigning CII missiles number 1 on the list. For the tasks 


allowing a choice of missiles, CII’s (the least capable missile) are given M? list position 


1, and CIIl’s position 2. The M? list assignment is shown in Table 3. 
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TASK GROUP M” List 
Tasks 1-13 
Tasks 14-26 
Tasks 27-39 
Tasks 40-52 
Tasks 53-65 
Tasks 66-78 
Tasks 79-91 
Tasks 92-104 
Table 3: M? List. Three task groups require CII missiles; all other task groups 
allow either CII or CII missiles. Where allowed, the less capable CII missiles are 
assigned M? list position one, and CII missiles position two. 





> >| > a 
RTP PMR TRIM ya 
= 


Scenario missile load-outs are identical for each ship of a particular class, and 
only two variants of TLAM are included. The cruisers are loaded with 32 TLAM’s as 
follows: 16 CII missiles, with one in Cell 1 of each of the 16 iagauiles. and 16 CIIl’s, one 
in Cell 5 of each module. ARLEIGH BURKE class destroyers are assigned 24 TLAM’s, 
12 CII’s, loaded in Cell 1 of each module, and 12 CII’s, loaded in Cell 5 of each module. | 
SPRUANCE class destroyers carry a full load of 61 TLAM’s. CII missiles are loaded 
into Cells 3,4,5, and 6; CIIl’s are waieiail to Cells 1,2,7, and 8. Figure 11 shows a 


representative load plan for one VLS launcher. 
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Sample CG Load-out 


s far [om 1 


Figure 11: Sample Cruiser Load Plan. All ships of a particular class are assigned 
identical load-outs. This figure shows an example of the load plan for one launcher 
of a TICONDEROGA class cruiser. Cruisers were loaded with a total of 16 CII 
missiles and 16 CIII missiles. 





B. RESULTS 

Table 4 gives a summary of scenario results. Results are provided for the 
weighted objective function, hierarchical restriction, and heuristic implementations where 
appropriate. The objective function values for each term are shown, as well as the total 
solution time. If an operator-selected option, such as primary missile spreading, or 
designating an expend platform, is not used for a scenario, then “N/A” is recorded for 


that objective function term. 


50 











ai ie pene tes Ae BY oe OES stage 


SCENARI 


Bn 9 ed RO ee 
HR | O| 0} NA| 34{ 3] NAT 13/8/81] 160] 26 | 
Li Heuristic | 0} 0] N/A] 34] 3] NAT 13/8 159[ 8 | 




























ST 8 OD 
4A BRO} NA} 10} 8.67] NAT N/A | 13/8/8 









Ne a _ 26/26/26 |__O| __ 2716 | 
2AHR | O| 1 NAT 425] 21 N/A| 26/20/16 
2A Heuristic | O] it N/A | 42.5 | 2 | N/A | 26/19/19 | 









Fe 1605" 
3Heuristic | 0| 0] NA] _155| 3] N/A] 52/32/32 “176. 







4uR_ | of it wat of 4] wal 65/52/40 | 
4 Heuristic | 0] NAL ol 4 | NA] 65/52/48 | 












SHR | oO] lo] NAT 114[ 4] AT 78/50/48 | 268 | 10976. 
SHeuristic | O| 0] NA| 114] 4] Nal 78/50/48 | 



















SAHR [| of of 65]  34[ 6] NA] 78/52/48 
Ce ee 













SBHR [| of 1[ NA[ 2147 4] Nal 78/62/48 | 
5B Heuristic | O} 1] N/A 214] 4] N/A 78/62/48 | 








a Od T 124/84/80 
6Heuristic | Of O] o69[ 24] 6] N/A| 126/8478|[ 204[ 98 


Table 4: Scenario Results. The first two base scenarios and their associated 





modifications have been solved using the weighted objective function method; all 


scenarios are solved with the hierarchical restriction and heuristic 


implementations. Objective function values are given for each sub-problem. For 
sub-problem 1 7, which minimizes the sum of M? list positions, results are given as 
the sum of M? list positions for primary, ready-spare, and back-up task parts. See 
Table 1 for the correspondence between objective function term and sub-problem 


number. 
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The weighted objective function method solution for sub-problem 7 is always 
poorer than either the hierarchical restriction or heuristic methods; this is an artifact of 
positive integrality gaps discussed in Chapter II. For the same reason, sub-problem 8, 
which maximizes residual salvo capability, is never considered with the weighted 
objective function method. A more detailed discussion of the positive integrality gaps 
demonstrated by this is presented later in this chapter. The weighted objective function 
method exhibits rapid growth of solution times with increasing problem size; in fact, this 
method did not return a solution for any scenario larger than Scenario 2. | 

We conclude that the weighted objective function is intractable for this 
application. 

A comparison of the solution times and objective Racten values between the 
hierarchical restriction and heuristic implementations provides insight into the relative 
merits of the heuristic. The scenarios presented here yield solution times for the heuristic 
that are always less than those obtained with the hierarchical restriction method. The 
absolute sence in solution times becomes more pronounced as the scenarios become 
more difficult to solve; this is discussed later in this chapter. The heuristic returns the 
same objective function value as the hierarchical restriction method for the first four sub- 
problems of every scenario tested. A discrepancy between the heuristic and hierarchical 
restriction objective function values is not realized until sub-problem 5, which is not 
optimized with the heuristic implementation, but rather evaluated using the fixed missile 
allocation from sub-problem 4. Of the 11 scenarios tested, only one, Scenario 5A, shows 


a difference between the objective function values. The hierarchical restriction yields a 
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superior solution for the only scenario that included back-up spreading. In this case, the 
heuristic also fixes the allocations to compute the objective function value. Only three 
scenarios (2A, 4, and 6) show a difference between the two methods for sub-problem 7, 
which minimizes the sum of M? list positions. These are particularly restrictive cases, in 
which additional considerations make it more challenging for the solver to find a solution 
within the specified optimality criteria; Scenario 2A divides four available firing 
platforms between two launch areas, Scenario 4 contains the most tasks of the scenarios 
with only four firing platforms, and Scenario 6 requires allocations made in Scenario 5 to 
be fixed while new tasking is added. The solution returned by the heuristic for these more 
difficult cases illustrates that further testing is required to ctennine the most robust way 
of considering M> list position. Except in one case, the heuristic method always returns a 
poorer solution to sub-problem 8, which maximizes residual salvo capability; however, 
the worst case is within 5% of the hierarchical restriction solution. 

The number of variables for each scenario is given in Table 5 for each sub- - 
_ problem for the hierarchical restriction and the heuristic. The weighted objective function 
method is solved as a monolithic problem; the number of variables for the entire model is 
given in the first column, where applicable. The number of discrete variables is reported 
rather than the total number of variables; this prevents variables whose values have been 
fixed in a previous sub-problem from being included in the number shown for the 
heuristic, and is done to allow a meaningful comparison between the relative problem 
sizes of the hierarchical restriction and heuristic implementations. The total number of 


_ variables in each scenario 1s within 10% of the number of discrete variables. Table 5 
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highlights the reduction in problem size that the heuristic realizes by evaluating a fixed 


allocation rather than optimizing sub-problems 5, 6, and 8. Additionally, Scenarios 1, 1A, 
1B, 2, 5, and 6 demonstrate the size of problem that can be omitted by not considering 


sub-problem 2, where possible. 
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Table 5: Number of Discrete Variables by Sub-Problem. For the orn 
implemented with weighted objective functions, the number of discrete variables 
for the entire model is shown in the first column. The number of discrete variables 
is Shown rather than the total number of variables. This prevents counting 
variables whose values have been fixed under the heuristic method. The total 
number of variables for each scenario is within 10% of the number of discrete 
variables for each method. 






Table 6 gives the number of constraints for each sub-problem for the hierarchical 


restriction and the heuristic; the total number of constraints for the weighted objective 


55 








function models is shown in the first column, where appropriate. Note there is no 
appreciable difference in the reported number of constraints between the hierarchical 
restriction and heuristic models; however, in actuality, with the heuristic method many of 
these equations are reduced to expressing that a 7 is equal to a constant, which 
implies that the number of coefficients in the constraint matrix and/or the total number of 


“real” constraints is greatly reduced for the heuristic implementation. 
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Table 6: Number of Constraints by Sub-Problem. For the scenarios implemented 
with weighted objective functions, the number of constraints for the entire model 
is shown in the first column. There is no appreciable difference in the number of 
equations reported for the hierarchical restriction and heuristic methods; 
however, in actuality, many of these constraints may be reduced to stating that a 
constant is equal to a constant. This leads to substantial reduction in the number 
of ‘real’ constraints that must be considered with the heuristic method. 






Table 7 shows the relative integrality gap achieved for each sub-problem. A 
maximum allowed relative gap of 5% was used with the weighted objective function. 


This results in absolute gaps on the order of 5.6x10''. As discussed in Chapter II, the 


a 





effect of this gap is seen in the sub-optimal values for variables contained in the lower 
priority objective function terms. The model selects missiles that result in a sub-optimal 
solution for considerations such as the sum of M’ positions and residual salvo capability; 
this is further indication that a weighted objective function is not well suited for the 
predesignation application. For the hierarchical restriction method, we observe that, for 
some sub-problems, almost all of the solution time is spent trying to find a feasible 
solution. Once this initial solution is found, an optimal solution is returned in a matter of 
seconds. This encourages us to experiment with an allowed relative gap of 100%, with 
the intention of having the branch and bound process terminate as soon as a feasible 
solution is found. We find that, for most scenarios, this results in negligible returned 
integrality gaps; based on this we do not believe that solution times would be adversely 
affected by specifying a (reasonable) tolerance on the order of 5%. This has not been 
done because, in almost all cases, a tolerance of 100% resulted in returned gaps of less 
than 5%. (The exceptions were Scenario 5B, nipaetieas 2, and Scenario 6, sub-problem 
1; we set the tolerance at 10%.) An integrality gap greater than zero only occurs in some 
instances of sub-problem 4, which minimizes the difference between the number of 
missiles remaining on a firing platform and the mean number of missiles on all firing 
platforms. Even in this case, the largest gap we observe is 5%; this equates to a difference 
of one missile. Note that because the hierarchical restriction and heuristic both solve sub- 
problem 4 in its entirety, the integrality gaps and objective function values for this sub- 


problem are identical for both methods for all scenarios. 
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Table 7: Returned Relative Integrality Gap by Sub-Problem. All gaps are no more 
than 0.05 (i.e., 5%). 


Solution times for each sub-problem are reported in Table 8. It is apparent that the 
heuristic takes longer to solve sub-problem 1 than the hierarchical restriction; this is 
because the heuristic attempts to first meet all allocations without using any platforms 
with employment penalties. If the required tasking cannot be met with this restricted set 


of available platforms, the sub-problem is re-solved using all of the firing platforms. 
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Additionally, the heuristic will attempt to solve the sub-problem for separate task parts, 
fixing the allocations after each task part is solved for. This may not allow a feasible 
solution, in which case the sub-problem is re-solved for all task parts together. In 
practice, this can lead to multiple attempts to solve sub-problem 1: (i) by individual task 
parts without considering firing platforms with employment penalties, (11) for all task 
parts simultaneously, with the restricted set of platforms, (iii) by individual task parts 
with all firing platforms considered, and finally (iv) for all task parts simultaneously, 
considering all firing platforms. However, the heuristic does reduce the total solve times 
in the instances where it omits sub-problem 2. Scenarios that do not designate expend 
platforms permit the elimination of sub-problem 3, reducing solve time. The full 
hierarchical restriction for sub-problem 4 is solved by the heuristic, and the solution times 
for the two methods are, in general, comparable. Evaluating sub-problems 5, 6, and 8 
using a fixed allocation provides a significant reduction in solution time. The heuristic 
attempts to solve sub-problem 7 for individual task parts, as it does for sub-problem 1. 
For the scenarios here, the results are inconclusive. if instances where a feasible solution 
can be found by solving for primary, ready-spare, and back-up task parts separately and 
fixing allocations after each task part, the heuristic takes less time. In cases where the 
heuristic cannot find a feasible solution, the entire sub-problem is re-solved for all task 
parts together, and the solution time for sub-problem 7 is greater than for the hierarchical 
_ restriction implementation. 

Unacceptably high solution times (e.g., about 44,000 seconds for one sub-problem 


instance) encouraged us to explore alternate executions of the CPLEX branch-and-bound 
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algorithm. To this end, care has been taken to achieve the best possible performance from 
CPLEX for all executions of the monolithic, hierarchical restriction, and heuristic 
methods, Extensive testing has determined the CPLEX parameter settings that result in 
the fastest execution of the branch-and-bound algorithm for each set of problem 
instances. For all monolithic problems, and for (both the hierarchical restriction and the 
heuristic implementation of) sub-problems 1 through 6, the settings that yield the best 
average performance are: i) select the branching variable based on “pseudo reduced 
costs’, (i1) solve the initial LP relaxation with the primal Simplex method, (iii) solve the 
sub-problems at each node with the primal Simplex method, (iv) choose the next node in 
the branch-and-bound tree for processing based on the best estimate value (given an 
integer completion) of the integer objective function, and (v) turn off the intrinsic CPLEX 
heuristic to realize feasible integer completions. Scenario 5B, because of the three ship . 
employment penalties, is more difficult to solve than the corresponding sub-problem for 
the other scenarios. In this instance, the best settings are: (i) branch upward first at each 
node, (ii) select the branch variable based on solving a sample of sub-problems to 
determine which branching variable might yield the best results, and (iii) use probing to 
tighten bounds at the root node (e.g., by fixing certain variable values). The best settings 
for sub-problem 7 include all of the above settings, but rather than selecting the 
branching variable based on “pseudo reduced costs”, the selection is based on “pseudo 
costs”. The best settings for sub-problem 8 include all parameter settings used for the 
first six sub-problems, and the branching and probing strategies used to solve sub- 


problem 7. [CPLEX, 1995] 
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Table 8: Solution Time by Sub-Problem. The ar coer in om solve 

‘times for sub-problem 1; the benefit of the heuristic is realized where the method 
of solving sub-problem 1 allows sub-problem 2 to be omitted. A further reduction 
is realized by the heuristic in cases that do not require an expend platform, 
allowing sub-problem 3 to be excluded. Evaluating sub-problems 5, 6, and 8 using 
a fixed allocation produces the greatest difference in total solution times between 
the hierarchical restriction and the heuristic. Attempting to solve sub-problem 7 
for individual task parts with the heuristic produces inconclusive results. In cases 
where a feasible solution can be found by solving for primary, ready-spare, and 
back-up parts in succession, the heuristic is faster. If a feasible solution cannot be 
found using this method, then the sub-problem is re-solved for all task parts 
together, and the hierarchical restriction is faster than the heuristic. 
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IV. CONCLUSIONS 


A. GENERAL 

Mixed integer programming does not lend itself readily to Tactical Decision Aids 
in general, and the TSC predesignation application in particular. There is no guarantee 
that a mixed integer program that has performed well for certain problem instances will 
perform adequately for alternate scenarios. The time required to reach a solution can vary 
widely as input data changes, and integer programming solvers (with tailored parameter 
settings) that appear adequate for one problem instance may prove disastrous when 
applied to another. These concerns complicate adoption of mixed integer programming as 
the underpinning of a tool that is intended for use by an operator in the fleet, who requires 
_ auseable solution in a reasonable amount of time to execute National Command 
Authority tasking. Through our implementation methods, we have attempted to provide 
the best possible structure for the use of an imperfect tool; it must be noted that each 
method has significant shortcomings that must be considered before the form of the TSC 
predesignation model is finalized. 


B. | WEIGHTED OBJECTIVE FUNCTION IMPLEMENTATION 


The weighted objective function is unsuitable for the TSC application because of 
the complication induced in integrality gaps and numerical stability. Poor solutions and 
rapid increase of solution times with problem size are unacceptable. The weights used for 
" our scenarios were computed using assumptions about the maximum number of firing 
platforms, weapons, and tasks. To use the weighted objective function in a real-world 


operational context, the weights would have to be modified to reflect the actual force 
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composition and level of tasking for each scenario; the validity of the model would 
depend on the extent to which the operator could perform these modifications. 
C. HIERARCHICAL RESTRICTION AND HEURISTIC 

Our hierarchical restriction accurately captures the desired priorities, and 
produces the best obtainable (although not necessarily optimal) solution to the 
predesignation problem. However, this method requires long solution times even for 
moderate-size scenarios. As strike weapons are developed that have preparation times 
substantially shorter than TLAM, the on-scene decision maker will need a tool that is 
capable of rapidly returning a solution. Even for the TLAM-only scenarios presented 
here, the hierarchical restriction does not appear to be sufficiently responsive. 

The heuristic returns solutions of reasonable quality in : significantly reduced 
amount of time for all the scenarios aed This method also recognizes all priority 
levels. However, the heuristic presented here is a first step in an evolution to a product 
robust enough for operational use. There are technical and doctrinal challenges that must 
be resolved before a heuristic is selected for use in the TSC predesignation system. In 
worst-case scenarios, the heuristic may need to solve the full hierarchical restriction for a 
sub-problem after several attempts to solve a restricted version of that sub-problem. This 
can lead to solution times that are longer for the heuristic than for the seminal 
hierarchical restriction. Additionally, the manner in which the heuristic evaluates fixed 
allocations for certain objective functions, rather than optimizing, needs to be clearly 
presented to the seat users of this product. The ordering a the priorities presented in 


this thesis has been developed jointly with NSWCDD [Fennemore, 1998c], and 








informally reviewed by members of Tactical Training Group Atlantic [Williams, 1998]; 
however, these orderings do not reflect official guidance, and will require doctrinal 


approval. 


D. RECOMMENDATIONS 


We recommend that NSWCDD develop a non-integer programming-based 
heuristic for TSC predesignation. However, a signal disadvantage of such a heuristic is 
that it is seldom possible to objectively assess the quality of its solutions. The 
mathematical models and formal methods introduced here may be useful in 
independently assessing the quality of the non-integer programming-based heuristic 
solution. Due to the limitations of mixed integer programming and of the heuristic 
presented here, we further recommend that at an early stage of development, the heuristic 
be presented to the strike warfare community for review. The goal of this review by 
decision makers is to determine its suitability for tactical systems, and to determine the 
correct ordering of priorities. We recommend that NSWCDD advocate the acceptance of 
some heuristic by the Naval Doctrine Command, Fleet Commanders, and Theater 
CINC’s, and its inclusion in the Naval Warfare Publications and Concepts of Operations 
that govern the conduct of strike warfare. It is crucial that strike planning be standardized 


in published doctrine so that it may be carried out with the effectiveness that it deserves. 
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